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Abstract 

Reaching a theoretical accuracy in the prediction of the lightest MSSM Higgs-boson 
mass, My, at the level of the current experimental precision requires the inclusion 
of momentum-dependent contributions at the two-loop level. Recently two groups 
presented the two-loop QCD momentum-dependent corrections to My mm , using 
a hybrid on-shell-DR scheme, with apparently different results. We show that the 
differences can be traced back to a different renormalization of the top-quark mass, 
and that the claim in Ref. [2] of an inconsistency in Ref. [lj is incorrect. We furthermore 
compare consistently the results for My obtained with the top-quark mass renormalized 
on-shell and DR. The latter calculation has been added to the FeynHiggs package and 
can be used to estimate missing higher-order corrections beyond the two-loop level. 
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1 Introduction 


The particle discovered in the Higgs-boson searches by ATLAS |3j and CMS [3] at CERN 
shows, within experimental and theoretical uncertainties, properties compatible with the 
Higgs boson of the Standard Model (SM) [5]- 7]- It can also be interpreted as the Higgs boson 
of extended models, however, where the lightest Higgs boson of the Minimal Supersymmetric 
Standard Model (MSSM) [8] is a prime candidate. 

The Higgs sector of the MSSM with two scalar doublets accommodates five physical 
Higgs bosons. In lowest order these are the light and heavy CP-even h and H, the CP-odd 
A, and the charged Higgs bosons H ± . At tree level, the Higgs sector can be parameterized in 
terms of the gauge couplings, the mass of the CP-odd Higgs boson, Ma, and tan/3 = V 2 /V 1 , 
the ratio of the two vacuum expectation values; all other masses and mixing angles follow 
as predictions. 

Higher-order contributions can give large corrections to the tree-level relations mum, 
and in particular to the mass of the lightest Higgs boson, M h . For the MSSM0 with real 
parameters the status of higher-order corrections to the masses and mixing angles in the neu¬ 
tral Higgs sector is quite advanced, see Refs. [T61 - 1T9] for the calculations of the full one-loop 
level. At the two-loop level [ 1511201 - 1341] in particular the 0(a t a s ) and (9(ctf) contributions 
(a t = h|/(47r), h t being the top-quark Yukawa coupling) to the self-energies - evaluated 
in the Feynman-diagrammatic (FD) as well as in the effective potential (EP) method - as 
well as the 0(aba s ), 0(a t ab ) and 0(al) contributions - evaluated in the EP approach - 
are known for vanishing external momenta. An evaluation of the momentum dependence 
at the two-loop level in a pure DR calculation was presented in Ref. [35]. The latest status 
of the momentum-dependent two-loop corrections will be discussed below. A (nearly) full 
two-loop EP calculation, including even the leading three-loop corrections, has also been 
published [36] . Within the EP method all contributions are evaluated at zero external mo¬ 
mentum, however, in contrast to the FD method which in principle allows for non-vanishing 
external momenta. Furthermore, the calculation presented in Ref. [3B] is not publicly avail¬ 
able as a computer code for Higgs-boson mass calculations. Subsequently, another lead¬ 
ing three-loop calculation of 0(a t a^), depending on the various SUSY mass hierarchies, 
was completed [37], resulting in the code H3m which adds the three-loop corrections to the 
FeynHiggs [TT1 !2T113&H40] result. Most recently, a combination of the full one-loop result, 
supplemented with leading and subleading two-loop corrections evaluated in the FD/EP 
method and a resummation of the leading and subleading logarithmic corrections from the 
scalar-top sector has been published [30] in the latest version of the code FeynHiggs. 

The measured mass value of the observed Higgs-boson is currently known to about 
250 MeV accuracy [5], reaching the level of a precision observable. At a future linear collider 
(ILC), the precise determination of the light Higgs-boson properties and/or heavier MSSM 
Higgs-bosons within the kinematic reach will be possible [4Tj. In particular, a mass mea¬ 
surement of the light Higgs-boson with an accuracy below ~ 0.05 GeV is anticipated [42], 

In Ref. [39j the remaining theoretical uncertainty in the calculation of Mh, from unknown 
higher-order corrections, was estimated to be up to 3 GeV, depending on the parameter 

*We concentrate here on the case with real parameters. For the case of complex parameters 
see Refs. [TTHTRj and references therein. 
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region; see also Refs. [3UJI33] for updated results. As the accuracy of the M h prediction 
should at least match the one of the experimental result, higher-order corrections which do 
not dominate the size of the ffiggs-boson mass values have to be included in the ffiggs-boson 
mass predictions. 

To better control the size of momentum-dependent contributions, we recently presented 
the calculation of the 0(p 2 a t a s ) corrections to Mh (the leading momentum-dependent two- 
loop QCD corrections). The calculation was performed in a hybrid on-shell/DR scheme [Ij 
at the two-loop level, where iff 4 and the tadpoles are renormalized on-shell (OS), whereas 
the ffiggs-boson fields and tan/3 are renormalized DR. At the one-loop level the top/stop 
parameters are renormalized OS H Subsequently, in Ref. |2j this calculation was repeated with 
a different result (also, a calculation in a pure DR scheme as well as the two-loop corrections 
of 0(aa s ) were presented). Within Ref. [2] the discrepancy between Refs. [1] and [2] was 
explained by an inconsistency in the renormalization scheme used for the ffiggs-boson held 
renormalization in Ref. [IJ. 

In this paper we demonstrate that this claim is incorrect. The renormalization scheme 
for the ffiggs-boson fields used in Ref. [I] is (up to corrections beyond the two-loop level) 
identical to the one employed in Ref. |2J. We clarify that the differences between the two 
results originates in a difference of the top-quark-mass renormalization scheme. While in 
Ref. [I] a full OS renormalization was used, in Ref. [2] the contributions to the top-quark 
self-energy of 0(e) (with 4 — D = 2e, D being the space-time dimension) were neglected, 
leading to the observed numerical differences. We also demonstrate how this difference in the 
treatment of the contributions from the top quark mass can be linked to a difference in the 
two-loop held renormalization constant and explain why this difference should be regarded 
as a theoretical uncertainty at the two-loop level, which would be hxed only at three loop 
order. 

We further present a consistent calculation of the 0(p 2 a t a s ) corrections to Mh in a 
scheme where the top quark is renormalized DR, whereas the scalar tops continue to be 
renormalized OS. This new scheme is available from FeynHiggs version 2.11.1 on, allowing 
for an improved estimate of (some) unknown higher-order corrections beyond the two-loop 
level originating from the top/stop sector. 

The paper is organized as follows. An overview of the relevant sectors and the renormal¬ 
ization employed in our calculation is given in Sect. [2] I 11 Sect. [3] we compare analytically 
and numerically the results of Refs. [T] and [2J. Results obtained using the DR scheme for 
the top-quark mass are given in Sect. [U Our conclusions are given in Sect. El 

2 From a technical point of view we calculated the momentum-dependent two-loop self-energy diagrams 
numerically using the program SecDec |44H46| . 
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2 The relevant sectors and their renormalization 


2.1 The Higgs-boson sector of the MSSM 

The MSSM requires two scalar doublets, which are conventionally written in terms of their 
components as follows, 

%= ( n ) = {v 2 + ^W + iX° 2 )) ' 

The bilinear part of the Higgs potential leads to the tree-level mass matrix for the neutral 
CP-even Higgs-bosons, 


M 2 ,tree _ 
VW Higgs - 


m 


0 i 


m 


0102 


m 


0102 


m 


02 


M\ sin 2 /? + Mf cos 2 /? - (M\ + Mf) sin /? cos (3 
— ( M\ + Mf) sin (3 cos /3 M\ cos 2 (3 + Mf sin 2 /? 


( 1 ) 


in the (</>!, 0 2 ) basis, expressed in terms of the Z boson mass, M z , Ma and the angle f3. 
Diagonalization via the angle a yields the tree-level masses m^t re e and m^ tree . Below we 
also use Mw, denoting the W boson mass and s w , the sine of the weak mixing angle, 

Sw = = y/1 


The higher-order-corrected CP-even Higgs-boson masses in the MSSM are obtained from 
the corresponding propagators dressed by their self-energies. The calculation of these and 
their renormalization is performed in the (</>!, <f 2 ) basis, which has the advantage that the 
mixing angle a does not appear and expressions are in general simpler. The inverse propa¬ 
gator matrix in the {(p \, 0 2 ) basis is given by 


(Affiggs) - -i 


P 2 -™0! +%(P 2 ) 

“ m 0i02 +^0i0 2 (P 2 ) 


_m 0102 + ^0102(P 2 ) ] 

P 2 -< + £ 02 (p 2 V 


( 2 ) 


where £(p 2 ) denote the renormalized Higgs-boson self-energies, p being the external mo¬ 
mentum. The renormalized self-energies can be expressed through the unrenormalized self¬ 
energies, F(p 2 ), and counterterms involving renormalization constants dm 2 and SZ from 
parameter and field renormalization. With the self-energies expanded up to two-loop order, 
t = S (1) + E (2) , one has for the CP-even part at the z-loop level (/ = 1, 2), 


£/,V) 

4'Pp 2 ) 

s£(p 2 ) 


£j,V) + iz® (P 2 , 

s 2*,(P 2 ) - sz tL m l«h ~ Sm 7L • 

S2(P 2 ) + SZ‘3 (P 2 - "4> - Snita ■ 


(3a) 

(3b) 

(3c) 


At the two-loop level the expressions in Eqs. (J2D do not contain contributions of the type 
(1-loop) x (1-loop); such terms do not appear at 0(a t a s ) and hence can be omitted in the 
context of this paper. For the general expressions see Ref. [IB]. 
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Beyond the one-loop level, unrenormalized self-energies contain sub-loop renormaliza¬ 
tions. At the two-loop level, these are one-loop diagrams with counterterm insertions at the 
one-loop level. 


2.2 Renormalization 


The following section summarizes the renormalization worked out in Ref. [T] , based on 
Ref. [21]. The held renormalization is carried out by assigning one renormalization con¬ 
stant to each doublet, 


ni^(i + isz Hl )n u U2 ^{i + \5z H2 )U2, (4) 

which can be expanded to one- and two-loop order according to 

SZ n , = 6Z%1 + 6Z%1 , SZ H , = sz£l + 6Z™ . (5) 

The held renormalization constants appearing in ([3j) are then given by 




- A zh) 
~ dZj Hi J 




- A 7® 
~ d/j H 2 ’ 




= i(6Z . (i) 


Hi 


+ 5Z, 


Hl> 


( 6 ) 


The mass counterterms 5m 2 ^ in Eq. (J3J) are derived from the Higgs potential, including the 
tadpoles, by the following parameter renormalization, 

M\ -)• M\ + 5M 2 a (1) + 5M 2 } 2) , T x -)■ Ti + + 5t[ 2) , (7) 

M| -> M| + 5M 2(1) + <5M 2(2) , r 2 -> r 2 + hT 2 (1) + hT 2 (2) , 

tan/3 —> tan/3 (l + /) tan/3*- 1 ) + /)tan/3^) . 


The parameters 7\ and T 2 are the terms linear in c/q and </> 2 in the Higgs potential. The 
renormalization of the Z-mass Mz does not contribute to the 0{a s at) corrections we are 
pursuing here; it is listed for completeness only. 

The basic renormalization constants for parameters and helds have to be hxed by renor¬ 
malization conditions according to a renormalization scheme. Here we choose the on-shell 
scheme for the parameters and the DR scheme for held renormalization and give the expres¬ 
sions for the two-loop part. This is consistent with the renormalization scheme used at the 
one-loop level. 

The tadpole coefficients are chosen to vanish at all orders; hence their two-loop counter¬ 
terms follow from 


T$ + 5T$ = 0 , 


i.e. 


§t (2) = _ t (2). 


( 2) = 


( 2 ) 


( 8 ) 


where T^ 2 \ are obtained from the two-loop tadpole diagrams. The two-loop renormal¬ 
ization constant of the A-boson mass reads 


SM 2 } 2) = ReE&M) 


( 2 ) 


(9) 


in terms of the A-boson unrenormalized self-energy Haa- The appearance of a non-zero mo¬ 
mentum in the self-energy goes beyond the 0(a t a s ) corrections evaluated in Refs. [2Q ].!2Tl[26] . 
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For the renormalization constants 5Z- Hl , 5Z^ 2 and 5 tan / 3 several choices are possible, see 
the discussion in m- As shown there, the most convenient choice is a DR renormalization 
of 5 tan/3, 5Z ni and 5Z n 2 , which at the two-loop level reads 


SZ% = 6Z+ {2) 


Re S 


'( 2 ) 

<Pl 


div 


\p 2 =0 


(10a) 


r 7(2) _ r ^DR(2) 

6Zj H 2 ~ 6/j H2 


Re S 


'( 2 ) 

</>2 


div 


b 2 =0 


(10b) 


5 tan = 5 tan /3 DR ( 2) = ^^SZ^ — SZ^j . (10c) 


The term in Eq. (11 Pel) is in general not the proper expression beyond one-loop order even in 
the DR scheme. For our approximation, however, with only the top Yukawa coupling at the 
two-loop level, it is the correct DR form [48]. 

The two-loop mass counterterms in the renormalized self-energies (J3]) are now expressed 
in terms of the two-loop parameter renormalization constants, determined above, as follows, 


5m 2 } 2) = 5M ^ z; cos 2 (3 + SM2*’ sin 2 /? 


r 2(2) 


2 ( 2 ) 


-5TI 


(2) 


cos/?(l + sin 2 /?) + YT 2 


(2) 


2M w s w 2M\ys v 

+ 25 tan/?^ cos 2 /?sin 2 /? {M\ — M|), 


cos 2 /? sin /? 


SmTL = -( SM + 5M sin /? cos (5 


r 2(2) 


2 ( 2 ) \ 


<t> 1<I>2 


- 5Tl 


( 2 ) 


e 1 ,3fl XT^( 2 ) 


sin /? - ST? 


cos 3 /3 


2M w s w 2M\ys 

— 5 tan /?® cos (3 sin j3 cos 2(3 {M\ + M|), 

5mf 2] = SM^ sin 2 /3 + SM 2 ^ cos 2 /? 


ST 


( 2 ) 


—- sin 2 /? cos /? — 5T^ -— 


2Mws w ~ 2M\ys v 

— 25 tan /3^ cos 2 /3 sin 2 /? (M\ — M§). 


sin/?(l + cos 2 (3) 


(11a) 


(lib) 


( 11 c) 


The Z-mass counterterm is again kept for completeness; it does not contribute in the ap¬ 
proximation of 0(atQ! s ) considered here. 


2.3 Diagram evaluation 

Our calculation is performed in the Feynman-diagrammatic (FD) approach. To arrive at 
expressions for the unrenormalized self-energies and tadpoles at 0(ata s ), the evaluation of 
genuine two-loop diagrams and one-loop graphs with counterterm insertions is required. For 
the counterterm insertions, described in subsection 12.41 one-loop diagrams with external 
top quarks/squarks have to be evaluated as well, as displayed in Fig. [Q The calculation is 
performed in dimensional reduction [45], 

The complete set of contributing Feynman diagrams was generated with the program 
FeynArts [50] (using the model Ele including counterterms from Ref. [51]), tensor reduction 
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and the evaluation of traces was done with support from the programs FormCalc [[52] and 
TwoCalc [53], yielding algebraic expressions in terms of the scalar one-loop functions A 0 , 
Bo |54j . the massive vacuum two-loop functions [55], and two-loop integrals which depend 
on the external momentum. These integrals were evaluated with the program SecDec [441 [46] . 
where up to four different masses in 34 different mass configurations needed to be considered, 
with differences in the kinematic invariants of several orders of magnitude. 



Figure 1: Generic one-loop diagrams for subrenormalization counterterms for the top quark 
(upper row) and for the scalar tops (lower row) (i,j,k = 1,2). 


2.4 The scalar-top sector of the MSSM 

The bilinear part of the top-squark Lagrangian, 

Amass = - (4, 4) M t , (12) 

contains the stop-mass matrix 

M? l + m 2 + Mf cos 2/3 (Tf - Q t sl ) m t X t 

m t X t M? r + m 2 + Mf cos 2/3 Q t sl 

with 

X t = A t - fi cot (3 (14) 

where Q t and Tf denote the charge and isospin of the top quark, A t the trilinear coupling 
between the Higgs bosons and the scalar tops, and /i the Higgsino mass parameter. Below 
we use Msusy := Mf L = M^ R for our numerical evaluation. The analytical calculation was 
performed for arbitrary Mf L and Mf , however. Mf can be diagonalized with the help of 
a unitary transformation matrix Uf, parameterized by a mixing angle 0f, to provide the 
eigenvalues rri- and m 2 ^ as the squares of the two on-shell top-squark masses. 

For the evaluation of the 0(a t a s ) two-loop contributions to the self-energies and tadpoles 
of the Higgs sector, renormalization of the top/stop sector at 0(a s ) is required, giving rise 
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to the counterterms for sub-loop renormalization. We follow the renormalization at the one- 
loop level given in Refs. [23 ]15SH58] . where details can be found. In particular, in the context 
of this paper, an OS renormalization is performed for the top-quark mass as well as for the 
scalar-top masses. This is different from the approach pursued, for example, in Ref. [35] . 
where a DR renormalization was employed, or similarly in the pure DR renormalization 
presented in Ref. [2], Using the OS scheme allows us to consistently combine our new 
correction terms with the hitherto available self-energies included in FeynHiggs. 

Besides employing a pure OS renormalization for the top/stop masses in our calculation, 
we also obtain a result in which the top-quark mass is renormalized DR. This new top-quark 
mass renormalization is included as a new option in the code FeynHiggs. The comparison 
of the results using the DR and the OS renormalization allows to estimate (some) missing 
three-loop corrections in the top/stop sector. 

Finally, at 0(a t a s ), gluinos appear as virtual particles only at the two-loop level (hence, 
no renormalization for the gluinos is needed). The corresponding soft-breaking glnino mass 
parameter M 3 determines the gluino mass, rrig = M 3 . 

2.5 Evaluation and implementation in the program FeynHiggs 

The resulting new contributions to the neutral CP-even Higgs-boson self-energies, containing 
all momentum-dependent and additional constant terms, are assigned to the differences 

AE rf (p 2 ) = Ei 2 6 V)-E2>(0), ab={HH,hH,hh}. (15) 

These are the new terms evaluated in Ref. [I], included in FeynHiggs. Note the tilde (not 
hat) on lh 2 )(0) which signifies that not only the self-energies are evaluated at zero external 
momentum but also the corresponding counterterms, following Refs. [201 12T] , A finite shift 
AE(0) therefore remains in the limit p 2 —> 0 due to 5M^ 2) = ReS^(M^) being computed 
at p 2 = M\ in E®, but at p 2 = 0 in E* 77 for details see Eqs. ([9]) and (JTTJ) . For the sake of 
simplicity we will refer to these terms as 0(p 2 a t a s ) despite the M\ dependence. 

3 Discussion of renormalization schemes 

In this section we compare our results for the (D(p 2 a t a s ) contributions to the MSSM Higgs- 
boson self-energies, as given in Ref. [I] to the ones presented subsequently in Ref. [2]. We first 
show analytically the agreement in the Higgs field renormalization in the two calculations 
and discuss the differences in the m t renormalizations. We also present some numerical 
results in both schemes, demonstrating agreement with Ref. [2] once the 0{e) terms are 
dropped from the top-quark mass counterterm. 

Using an OS renormalization for the top-quark mass, the counterterm is determined from 
the components of the 0(a s ) top-quark self-energy (Fig. [T]) as follows, 

—— = \ Re { ( m t) + ( m t)] + [ s f L ("h 2 ) + )] } , (16) 

mt z 

where the top-quark self-energy is decomposed according to 

£ tip) =^-^t(P 2 )+ (P 2 ) + m t uj^ L (p 2 ) + m*u; + EfV) • ( 17 ) 
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Figure 2: One-loop subrenormalization diagram contributing to <5s 22 (p 2 ) and 5a{p 2 ), with 
the counterterm insertion denoted by a cross. The right diagram only contributes to ^s 22 (0) 
and 5a (0). 


with the projectors u± = |(1 ± 75 ). 

3.1 Analytical comparison 

In the 0(a t a s ) calculation of the Higgs-boson self-energies the renormalization of the top- 
quark mass at 0(a s ) is required. The contributing diagrams are shown in the top row of 
Fig. [Q The top-quark mass counterterm is inserted into the sub-loop renormalization of 
the two-loop contributions to the Higgs-boson self-energies, where two sample diagrams are 
shown in Fig. [21 The left diagram contributes to the momentum-dependent two-loop self¬ 
energies, while the right one contributes only to the momentum-independent part. Evaluat¬ 
ing the expression in Eq. (fl 6 l) in 4 — 2e dimensions yields the OS top-quark mass counterterm 
at the one-loop level, which can be written as a Laurent expansion in £, 

5mf s = - 5mf v + 5mf n + £ 5m e t + ... ; (18) 

higher powers in e, indicated by the ellipses, do not contribute at the two-loop level for 
e —y 0 after renormalization. Accordingly, the DR top-quark mass counterterm is given by 
the singular part of Eq. (|T 8 l) , 

5m P = i 5mf v . (19) 

For further use we define the quantity 

5m™ = i 5mf v + 5mf n . (20) 

At 0(a s ) the OS counterterm is given as 

5m° s 
m t 


8 


67 r 


An(m+) , 9 9 . 

2 o -4R 0 (m 2 ,0,m 2 ) 


mi 


n M m !) , M m i) , M m l) 

^ 9 T 9 + 9 


mi 


mi 


mi 


m~ + m\ — m~ — 4 sin Of cos Of mg m t 

m t 

m~ + m? — mf +4 sin Q~ t cos 6f m^mt 


mi 


Re[H 0 (m 2 ,m~,m~)] 


Re[R 0 (m f 2 ,m?,m|)] 


( 21 ) 













The one- and two-point, functions A 0 (m 2 ) and B 0 (p 2 , m 2 , m 2 ) are expanded in e as follows, 


A 0 (m) = -A^lm) + A^im) +e^4o( m ) > 

£ 

B 0 (p 2 , ml, m 2 2 ) = - B$ lv (p 2 , m\, m 2 ) + B^{p 2 , m\, m 2 2 ) + e B £ (p 2 , m\, m 2 ) . (22) 


Consequently, the term at O(e), 5m £ t /m t , is given by Eq. (ED), but taking only into account 
the pieces oc Aq, Bq. The special cases of A^m 2 ) and _Bg(m 2 ,0,m 2 ) are given by 


A £ 0 (m 2 ) =m 2 <J 1 - log(m 2 //i 2 ) + ^ log 2 (m 2 //i 2 ) + ^ f , 


1 


7T 


So(m , 0, m ) = 4 - 2 log(m //i ) + - log (m jn ) + — , 


(23) 


where the factor 47re lE is absorbed into the renormalization scale. The expression for Bq 
depending on three mass scales can be found e.g. in Ref. [59]. 

In our calculation in Ref. [1] we include terms up to 0(e), originating from the top-quark 
self-energy, in the top-mass counterterm[§, i.e. 

6mP = 5m° s . (24) 


The derivation in Ref. [2] proceeds differently. The renormalized Higgs-boson self-energies 
are first calculated in a pure DR scheme. This concerns the top mass, the scalar-top masses, 
the Higgs field renormalization, and tan/3. In this way it is ensured that in particular 
the Higgs fields are renormalized using DR, 5Z^. = where this quantity contains 

the contribution from the one- and two-loop level. Using this pure DR scheme a finite 
result is obtained in which all poles in 1/e and \je 2 cancel, such that the limit £ —y 0 
can be taken. Subsequently, the DR top-quark mass counterterm, <5mP R , is replaced by 
an on-shell counterterm, and the top-quark mass definition is changed accordingly. (The 
same procedure is applied for the scalar-top masses.) Since these finite expressions for the 
renormalized Higgs-boson self-energies do not contain any term of 0(1/e), the 5m £ t part of 
the OS top-quark mass counterterm does not contribute, i.e. 

5mp = 5mf IN . (25) 


The numerical results for the renormalized Higgs-boson self-energies obtained this way differ 
significantly from the ones obtained in Ref. [I], as pointed out in Ref. [2]. 

In the following we discuss the different Higgs-boson field renormalizations, where we use 
the notation of SZ^ for the field renormalization derived using Smf, with X = DR, FIN, 
OS. The field renormalization can be decomposed into one-loop, two-loop, .. .parts as 


5Z%f = + 


(26) 


3 Taking 0(e) terms into account in the expressions for on-shell counterterms beyond one loop is widely 
used in the literature, see e.g. Refs. [60H62] . 
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In Ref. PJ it was claimed that using an OS top-quark mass renormalization from the 
start results in a non-DR renormalization of SZ% 2 . While it is correct that an OS value for 
m t yields different results in the one- and two-loop part, 


sz: 


Sm° s (l) 

n 2 


± sz: 


5 m]P R (l) 


n 2 


SZ, 


5m t os (2) 
H 2 


^ SZ 


*n? R ( 2 ) 
U 2 ’ 


(27) 


the sum of the one- and two-loop parts are identical, independently of the choice of the 
top-quark mass renormalization (see e.g. Eqs. (3.60)-(3.62) in Ref. (63]), 


=) C 


Sm? s 


div 


x r^Sm™ c 7 < 5 m,p R 
= SZ H, = SZ M2 


(= . 


(28) 


provided that also in SZ all finite pieces are dropped, as done in Ref. [Ij. Differences 
between SZ^ 2 and SZ® 2 arise only at the three-loop level. Consequently, the claim in Ref. [2] 
that using Srn^ leads to an inconsistency in the Higgs held renormalization in Ref. [I] is not 
correct. The held renormalizations thus cannot be responsible for the observed differences 
between Refs. (I] and pj. 

More explicitly, the difference between the two calculations results from non-vanishing 
Srrf t terms in the renormalized Higgs-boson self-energies. Those terms naturally appear 
when performing a full expansion in the dimensional regulator e. The latter corresponds to 
choosing (5m° s (as done in Ref. m) instead of <5mf IN (as done in Ref. 0). 

In order to isolate the contributions coming from O(s) terms x 1/e poles we define the 
following quantities, where superscripts OS, FIN refer to the respective use of 5m t os , 5mf IN : 


St^)°s = 5T &™ + ^ ; ( 29 a ) 

sg o V)= £<rv) +^«(p 2 ) ■ < 29b ) 

S<T(P 2 ) = sS FIN (p 2 ) + S A (p 2 ) , (29c) 


where the last equation yields a shift for the H-boson mass counterterm in Eq. (I7|) , 

SMf )os = SM 2 } 2) fin + S a (M\) . (30) 


The 5-terms are defined as the finite contributions stemming from <5mf-dependent parts 
in the counterterms (see the left diagram in Fig. [2] for an example). The DR-renormalized 
quantities do not contain a Unite (fm^-dependent part by definition. Furthermore, since 0i 
has no coupling to the top quark, there are no terms proportional to Smf in , S^h, and 

(5T] (2) , and it is sufficient to consider <5 e 22 , <5^, and <5 t 2 only. While St 2 is p 2 - independent, we 
find 


<fe 2 2 ( P 2 ) 
Sa(p 2 ) 


3a t 2 Sm e t 

TTP — + <W 0 ), 

27T m t 

3 a t 9 9 Sm e . _ , . 

tt 1 P 2 cos 2 /3 —- + <5^(0) 

27t m t 


(31) 

(32) 
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Using Eqs. (|3]h (fTTj) we find that the following relations hold for the renormalized Higgs- 
boson self-energies: 


^£ 22 ( 0 ) 


- sin 2 /3 5^(0) - ——-cos 2 /3sin/3 5t 2 

2M w s„ 

sin j3 cos /3 5R0) H-cos 3 (3 5t 2 

2M w s w 

cos 2 ^ <5j(0) + sin/J(l + cos 2 /?) <5 t 2 


0 (for e//) , 

0 (forU/j. 

0 (for £&>). 


(33) 


This is in agreement with the observation that in the renormalized Higgs-boson self-energies 
at zero external momentum at 0(ata s ), the terms containing 5m e t drop out in the final 
(finite) result. Such a cancellation is to be expected as the same combination of one-loop 
self-energies that potentially contributes to this finite contribution also appears in the 0(1/e) 
term, where they must cancel. This argument in principle still holds when the momentum- 
dependent 0(a t a s ) corrections are calculated and all counterterms are evaluated with a full 
expansion in e. Since the counterterm 5a is evaluated at p 2 = M\ |, and the Higgs-boson 
fields are renormalized in the DR scheme, however, one finds, using Eqs. (J3j), (TTTT) for the 
three renormalized Higgs-boson self-energies, 

- sin 2 /3 ( 5 a (M\ ) - 5 a ( 0)) = (- cos 2 /3 sin 2 /3 Mj) — (for eJJ } ) , 

ztt m t 

+ sin [3 cos (3 ( 5 a (M 2 a ) - 5 A (0)) = (+cos 3 ^ sin/3 M\) ^ (for , 

(^£22 (P 2 ) - ^£ 22 ( 0 )) - cos 2 /3 (5 a (M\) - 5 a (0)) = ^(p 2 ~ cos 4 /3 M%) ^ (for E®), 

(34) 


i.e. the 5rrf t terms contribute in the newly evaluated 0(p 2 a t a s ) corrections. They are p 2 - 
independent in E^ 2) and E^ 2 , while they do depend on p 2 in S^ 2 j. 

The p 2 -dependent terms coming from the expansion of terms like (— p 2 )~ e multiplying 
a 1 /e 2 divergence must certainly cancel after inclusion of the counterterms, because non¬ 
local terms cannot appear in a renormalizable theory. However, the cancellation of the 
e-dependent terms stemming from the mass renormalization is not necessarily fulfilled once 
the two-loop amplitude carries full momentum dependence. Similarly, the truncation of 
the field renormalization to the divergent part cuts away terms involving 5mf, leading to 
further non-cancellations. The explicit DR renormalization of the Higgs-boson fields drops 
the corresponding finite contributions, such that no 5m/ n , 5m e t terms are taken into account. 
The different dependence on the external momentum and the DR prescription for the Higgs 
field renormalization leads to Eqs. (I34]h 

Equivalent momentum-dependent terms of 0(e) of the scalar-top mass counterterms, 
evaluated from the diagrams in the lower row of Fig. [U do not contribute. The diagrams with 
top-squark counterterm insertions are depicted in Fig. [3l The first diagram is momentum 
independent. In the second diagram, the corresponding loop integral is a massive scalar 
three-point function (Co) with only scalar particles running in the loop, and thus is UV 
finite. Consequently, the top-squark mass counterterm insertions of 0(e ) do not contribute. 
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Figure 3: One-loop subrenormalization diagrams containing top-squark loops with countert¬ 
erm insertions. 


In the third diagram the stop mass counterterm can enter via the (dependent) counterterm 
for A t [2TII5T] . This diagram does not possess a momentum-dependent divergence, however, 
and thus the O(e) term of the scalar top mass counterterm again does not contribute. 


3.2 Physics content and interpretation 

In the following we give another view on the finite 5m £ t term from the top mass renormaliza¬ 
tion and on the interpretation of the different results for the Higgs-boson masses with and 
without this term. 


In the approximation with p 2 = 0 for the two-loop self-energies, the results are the same 
for either dropping or including the 5m £ t term, provided that this is done everywhere in the 
contributions from the top-stop sector in the renormalized two-loop self-energies. 


As explained above, abandoning the p 2 = 0 approximation yields an additional 5m £ t in 
the p 2 -coefficient of the self-energy E^ (p 2 ) when the on-shell top-quark mass counterterm, 
see Eq. m, is used, as well as in the A-boson self-energy E aa(p 2 ) from which it induces an 
additive term ~ M 2 A bm £ t /mt to the mass counterterm 8M\. 

In the renormalized self-energy E^(p 2 ), Eq. (l3c]h this extra p 2 -dependent term survives 


when SZ^l is defined in the minimal way containing only the 1/e and 1/e 2 singular parts; 


however, it disappears in E^ 2 J(p 2 ) when the minimal SZ^ = bZ'^// 


( 2 ) _ £ ( 2 ) 


div 


is replaced by 




3 a t brrf t 
2i r m t ’ 


(35) 


which now accommodates also a finite part of two-loop order. 

This shift in bzff by a finite term has also an impact on the counterterm for tan /3 via 

<5tan/3 = ^SZy^. This has the consequence that the extra bm\ term in 5M\ drops out 
in the constant counterterms for the renormalized self-energies E^ (p 2 ) in Eq. (j3j) because 

of cancellations with the 5m £ t term in htan/3 and 5Z^ (this can be seen from the explicit 
expressions given in Eqs. (J6J) and (fTTf ). 

Accordingly, keeping or dropping the hnite 8m £ t part is thus equivalent to a hnite shift 
in the field-renormalization constant bZ-^ 2 at the two-loop level, which corresponds to a 
hnite shift in tan/3 as input quantity. Numerically, the shift in tan/3 is small, and cannot 
explain the differences in the M/, predictions from the two schemes. Hence, these differences 
originate from the different p 2 coefficients in E^(p 2 ). 
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The impact of a modification of the two-loop field-renormalization constant on the mass 
Mh can best be studied in terms of the self-energy E^ in the h. H basis, which is composed 
of the in the following way, 

Y>hh = cos 2 ct E^ + sin 2 Q! E^ — 2 sinct cosct E^^, (36) 

where only E ^ 2 contains the p 2 -dependent 5m\ contribution. In order to simplify the discus¬ 
sion and to point to the main features, we assume sufficiently large values of tan /3 that we 
can write E hh — and h, H mixing effects play only a marginal role (both simplifications 

apply to the numerical discussions in the subsequent section). Moreover, to simplify the 
notation, we drop the indices and define 

E hh = E, E hh = E, SZ h h = SZ , (37) 

where 8Z hh = cos 2 a5Z n2 + sin 2 a5Z ni ~ SZ n2 . Starting from the tree-level mass m h and 
the renormalized h self-energy up to the two-loop level, 

E(p 2 ) = E(p 2 ) - Sml + 5Z(p 2 - ml) , (38) 

we obtain the higher-order corrected mass Mh from the pole of the propagator, i.e. 

Ml -m\ + t(Ml) = 0 . (39) 

The Taylor-expansion of the unrenormalized self-energy around p 2 = 0, 

E(p 2 )= E(0) + p 2 E'(0) + S(p 2 ), (40) 

yields the first two terms containing the singularities in 1/e and 1 /e 2 , and the residual fully 
finite and scheme-independent part denoted by E (p 2 ). With this expansion inserted into 
Eq. (l38li one obtains from the pole condition Eq. (l39li the relation 

(Ml - m 2 h ) [1 +SZ + E'(0)] + [£(0) - Sm\ + m\ E'(0)] + t(Ml) = 0, (41) 

where the expressions in the square brackets are each finite, irrespective of a possible finite 
term in the definition of 5Z. 


Taking into account that M/ differs from ml by a a higher-order shift, we can replace 


t(Ml) = t(ml) + (Ml - ml) t '(ml) + ■■■ 


(42) 


and obtain 

M l ~ m l = 


E(0) — Sml + m 2 h £'(0) + S(m^) 
1 + 5Z + E'(O) + t'(ml) 


(43) 


= - [m - Sml + ml £'(0) + E(m|)] lloop + 21oop 

+ [E(0) - Sml + ml E'(0) + E(m 2 J] llopp • [SZ + £'(0) + t' (ml)] llppp + • • ■ 

showing explicitly all terms up to two-loop order. It does not contain the two-loop part 
of the field-renormalization constant, which indeed would show up at the three-loop level. 
Hence, effects resulting from different conventions for hZl 21oop ) in the finite part have to be 
considered in the current situation as part of the theoretical uncertainty. 


13 



3.3 Numerical comparison 

In this section the renormalized momentum-dependent 0(p 2 ata s ) self-energy contributions 
A T,hh, A 'EhH, A Tjhh of Eq. (fl5|) and the mass shifts 

AMh = M h — M h!0 , A M h = M h — M H) o (44) 

are compared using either 5m° s or 5m™, as discussed above. M/^o and Mup denote the 
Higgs-boson mass predictions without the newly obtained (D(p 2 atOt s ) corrections. 

The results are obtained for two different scenarios. Scenario 1 is adopted from the m^ ax 
scenario described in Ref. [64]. We use the following parameters: 

m t = 173.2 GeV, M SUS y = 1 TeV, X t = 2M SUS y , 

rrig = 1500 GeV, /i = M 2 = 200 GeV. (45) 

Here M 2 denotes the SU( 2) soft SUSY-breaking parameter, where the U( 1) parameter is 
derived via the GUT relation Mi = (5/3) (s^/c/) M 2 . Scenario 2 is an updated version of 
the “light-stop scenario” of Refs. [641165] 

m t = 173.2 GeV, Msusy = 0.5 TeV, X t = 2 Msusy , 

m- g = 1500 GeV, p = M 2 = 400 GeV M x = 340 GeV, (46) 

leading to stop mass values of 

m h = 326.8 GeV, m h = 673.2 GeV . (47) 

A renormalization scale of p = mt is set in all numerical evaluations. 

Self-energies 

In Fig. [4]we present the results for the 5a (upper plot) and <5 e 22 (lower plot) contributions 
for tan/3 = 5(20) in red (blue) in Scenario 1, where 5a,5-£ 22 are defined in Eqs. (J29]). In the 
upper plot 5a(M\) (<5a( 0)) is shown as solid (dashed) line; correspondingly, in the lower plot 
^ 22 G° 2 ) (^ 22 ( 0 )) is depicted as solid (dashed) line. The contribution is seen to decrease 
quadratically with M 4 or p (:= ■s/p 2 ) when including the momentum-dependent terms, see 
Eq. (l34]h For 5a it is suppressed with tan 2 f3. For high values of Ma and low tan/3, the 5a 
contribution becomes sizable. Similarly, for large p the 5y, 22 term becomes sizable, showing 
the relevance of the 5m £ t contribution. 

The behavior of the real parts of the two-loop contributions to the self-energies AE a b is 
analyzed in Fig. [5] Solid lines show the result evaluated with <3m/ )S , as obtained in Ref. [T] 
(i.e. the new contribution added to the previous FeynHiggs result in Ref. [I], see Eq. (fl5]b . 
Dashed lines show the result evaluated with 5m™, as obtained in Ref. [2]. We show Ma = 
250 GeV and tan/3 = 5(20) as red (blue) lines. The difference between the 5mf IN and <5m/ ,s 
calculations for AS ^ )1 and AS^ 19 i 2 is p- independent, as discussed below Eq. (I34lh and the 
difference between the two schemes is numerically small. For AE^ 2 , on the other hand, the 
difference becomes large for large values of p. This self-energy contribution is mostly relevant 


14 




Figure 4: 5a{M\) and <5^(0) varying Ma shown in the upper plot, 22 (p 2 ) an d fe 22 (0) i n 

the lower plot, both within Scenario 1. 

for the light CP-even Higgs-boson, however, i.e. for p ~ Mh, and thus the relevant numerical 
difference remains relatively small (but non-zero) compared to the larger differences at large 
p. For completeness it should be mentioned that the imaginary part is not affected by 
the variation of the top-quark renormalization, as only the real parts of the counterterm 
insertions enter the calculation. 

Scenario 2 was omitted as the relevant aspects for the analysis of the self-energies using 
5m^ s vs. 5mf m have become sufficiently apparent within Scenario 1. 
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Mass shifts 

We now turn to the effects on the neutral CP-even Higgs-boson masses themselves. The 
numerical effects on the two-loop corrections to the Higgs-boson masses Mh t H are investigated 
by analyzing the mass shifts A Mh and A Mh of Eq. (j44|i . The results are shown for the two 
renormalization schemes for the top-quark mass, i.e. using or 5mf IN . The color coding 

is as in Fig. [5] The results for Scenario 1 are shown in Fig. [6] and are in agreement with 
Figs. 2 and 3 (left) in Ref. [2], i.e. we reproduce the results of Ref. m using Smf - [N . The 
results for Scenario 2 are shown in Fig. [TJ The results are again in agreement with Figs. 2 
and 3 (right) in Ref. [2j. This agreement confirms the use of 5mf 1K in Ref. [2] , in comparison 
with <5m° s used in the evaluation of our results. For the contribution to Mh , peaks can be 
observed at Ma = 277+ rny 2 , 2 m^, see also Ref. [I] and the discussion of Fig. [9] below. 

Since the the results using 5m® s and 5mf IN correspond to two different renormaliza¬ 
tion schemes, their difference should be regarded as an indication of missing higher-order 
momentum-dependent corrections. 
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Figure 5: AE^. in Scenario 1 (with Ma = 250 GeV) for ij = 11,12,22 in the upper, the 
middle and the lower plot, respectively. The solid (dashed) lines show the result obtained 
with (<5mf IN ); the red (blue) lines correspond to tan/3 = 5(20). 
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M a (GeV) 

Figure 6: Variation of the mass shifts A Mh, A Mh with the A-boson mass Ma within Sce¬ 
nario 1, for tan (3 = 5 (red) and tan/3 = 20 (blue) in- or excluding some 6 terms. The peak 
in A Mh originates from a threshold at 2 m*. 
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Figure 7: Variation of the mass shifts AM/,, AM# with the A-boson mass M 4 within Sce¬ 
nario 2, for tan (3 = 5 (red) and tan /? = 20 (blue) in- or excluding some 6 terms. The peaks 
in AM# originate from thresholds at 2m t , 2 + m^ 2 and 2 , where the threshold 

at 2 m t is suppressed by 1 /tan 2 /3. 
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4 Comparison with the m t DR renormalization 

Having examined the renormalization of the top-quark mass, we will now analyze the nu¬ 
merical differences between an mf R and an m° s calculation. This has been realized by 
employing a DR renormalization of the top-quark mass in all steps of the calculation. The 
top-squark masses are kept renormalized on-shell. This can be seen as an intermediate step 
towards a full DR analysis. 


4.1 Implementation in the program FeynHiggs 

In the DR scheme the top-quark mass parameter entering the calculation is the MSSM DR 
top-quark mass, which at one-loop order is related to the pole mass mt (given in the user 
input) in the following way, 


m 


DR 


(/■0 = m ■ 




(48) 


The term 5mf n can be obtained from Eq. (I18h . with the formal replacement a s —> a^ R (/i), 
yielding 


5mf n 

m t 


= a 


DR/ 


5 1 . , 2 /.^ , m 9 


(h) -IT + - log (m t /n ) + 


37T 

1 


7T 


3m? 7T 


(-1 + log (m?//i 2 )) 


+ - Jog«//i 2 )) + mf 2 (l - log(m| 2 /p 2 )) 


(m~ + m 2 — m\ — 2 mgm t sin(2 9 t )) Re[Z?o n (m 2 , m~, m 2 )] 


+ (m? + m\ - m\ + 2 m- g m t sin(20 t )) Rm? m| )] 


(49) 


At zeroth order, a° R (/i) = a^ s (/i). 

As on-shell renormalized quantities the stop masses mq and mi should have fixed values, 
independently of the renormalization chosen for the top-quark mass. We compensate for the 
changes induced by 5mf n in the stop mass matrix, Eq. f[T3]) . by shifting the SUSY-breaking 
parameters as follows, 


+ K os ) 2 - ( m ? R ) 2 , 
Ml ^ Mg =M? r + (m° s ) 2 - (m? 5 ) 2 , 


A t —» A' t 


m 


os 


m p R 



h 

tan /3 


(50a) 

(50b) 

(50c) 


(Except for A t , which actually appears in the Feynman rules, FeynHiggs only pretends to 
perform these shifts but computes the sfermion masses using 

This procedure is available in FeynHiggs from version 2.11.1 on and is activated by 
setting the new value 2 for the runningMT flag. The comparison of the results with DR and 
with OS renormalization admits an improved estimate of (some) of the missing three-loop 
corrections in the top/stop sector. 
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4.2 Numerical analysis 

In the following plots we show the difference 

Ail^ := ih4(m° s ) - </> = h,H, (51) 

between evaluated in the OS scheme, i.e. using m° s ( not mf IN ), and in the DR scheme, 
i.e. using 

Dependence on Ma 

In the upper half of Fig. (8J A Mh is plotted in Scenario 1 as a function of Ma with tan j3 = 
5(20) in red (blue). The solid (dashed) lines show the difference evaluated at the full one- 
loop level (including the 0(a t a s ) corrections). The dotted lines include the newly calculated 
0(p 2 a t a s ) corrections. For Ma > 200 GeV one observes large differences of 0( 10 GeV) at 
the one-loop level, indicating the size of missing higher-order corrections from the top/stop 
sector beyond one-loop. This difference is strongly reduced at the two-loop level, to about 
~ 3 GeV, now corresponding to missing higher orders beyond two-loop from the top/stop 
sector. The dotted lines are barely visible below the dashed lines, indicating the relatively 
small effect of the 0(p 2 a t a s ) corrections as derived in Ref. [Tj. 

The lower plot of Fig. [8]shows the corresponding results for A Mh with the same color/line 
coding. Here large effects are only visible for low Ma , where the higher-order corrections to 
Mh are sizable (and the light Higgs-boson receives only very small higher-order corrections). 
In this part of the parameter space the same reduction of A Mh going from one-loop to 
two-loop can be observed. 

The behavior is similar for Scenario 2, shown in Fig. [9] (with the same line/color coding 
as in Fig. [8]), only the size of the difference A Mh is ~ 20% smaller at the one-loop, and 
~ 50% smaller at the two-loop level compared to Scenario 1. The same peak structure due 
to thresholds as in Fig. [7] is visible. 
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M a (GeV) 



M a (GeV) 

Figure 8: A = M^(mf s ) — M^(mf R ) for 4> — h (upper plot) and (j) — H (lower plot). 
The difference is shown as solid (dashed/dotted) line at the one-loop ( 0(a t a s )/0(p 2 a t a s )) 
level as a function of Ma for tan/3 = 5(20) in red (blue) within Scenario 1. 
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Figure 9: A— M^(rn[ m ) for (p — h (upper plot) and </> = H (lower plot) as 
a function of Ma within Scenario 2, with the same line/color coding as in Fig. [8] The peak 
in the lower plot originates from a threshold at 2 . The threshold at 2 m t is suppressed 

by 1/ tan 2 j3. 
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Dependence on rrig 

In Figs, m EH we analyze Aas a function of rrig in Scenario 1 and 2, respectively. We 
fix Ma = 250 GeV and use the same line/color coding as in Fig. |8j Due to the choice of an 
MSSM DR top-quark mass definition, mf R varies with rrig already at the one-loop level. 

In the upper plots we show the light CV-eve n Higgs-boson case, where it can be observed 
that the scheme dependence is strongly reduced at the two-loop level. It reaches 2 — 3 GeV 
in Scenario 1 and ~ 1 GeV in Scenario 2, largely independently of tan f3. At the one-loop 
level the scheme dependence grows with rrig, whereas the dependence is much milder at 
the two-loop level. The effects of the 0(p 2 ata s ) corrections become visible at larger rrig, in 
agreement with Ref. |T| . 

The heavy CV-eve n Higgs-boson case is shown in the lower plots. At small tan /3 scheme 
differences of 0(600 MeV(150 MeV)) can be observed at the one-(two-)loop level. For large 
tan (3 the differences always stay below 0(50 MeV), in agreement with Fig. [3 The depen¬ 
dence on rrig is similar as for the light Higgs-boson, but again somewhat weaker. 

Dependence on X t 

Finally, in Figs. Q21 [13 we analyze Aas a function of X t = X t os in Scenario 1 and 2, 
respectively. We again fix Ma = 250 GeV and use the same line/color coding as in Fig. [3 

In the upper plots we show the light CV-eve n Higgs-boson case. As before the scheme 
dependence is strongly reduced when going from the one-loop to the two-loop case. In 
general a smaller scheme dependence is found from small X t , while it increases for larger 
\X t \ values, in agreement with Ref. [66]• For most parts of the parameter space, when the 
two-loop corrections are included, it is found to be below ~ 3 GeV. The contribution of 
0(p 2 a t a s ) remains small for all X t values. 

In the heavy CV-eve n Higgs-boson case, shown in the lower plots, the dependence of the 
size of the effects is slightly more involved, though the general picture of a strongly reduced 
scheme dependence can be observed here, too. In both scenarios, for large negative X t 
and tan/3 = 5 the 0(p 2 a t a s ) contributions can become sizable with respect to the 0(a t a s ) 
corrections. 

In conclusion, the scheme dependence is found to be reduced substantially when going 
from the pure one-loop calculation to the two-loop 0(ata s ) corrections. This indicates 
that corrections at the three-loop level and beyond, stemming from the top/stop sector are 
expected at the order of the observed scheme dependence, i.e. at the level of ~ 3 GeV. This 
is in agreement with existing calculations beyond two-loop mm- A further reduction of 
the scheme dependence might be expected by adding the 0(a 2 ) contributions. The mf R 
value calculated at 0(a s + a t ) is substantially closer to m ( os , reducing already strongly the 
scheme dependence at the one-loop level. This extended analysis is beyond the scope of our 
paper, however. 
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Figure 10: Afor 0 = h (upper plot) and <fi = H (lower plot) as 
a function of rrig within Scenario 1, for M 4 = 250 GeV and with the same line/color coding 
as in Fig. 0 
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Figure 11: A— M^m^) for 0 = h (upper plot) and <fi = H (lower plot) as 
a function of ntg within Scenario 2, for M 4 = 250 GeV and with the same line/color coding 
as in Fig. 0 
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Figure 12: Afor 4> — h (upper plot) and (f> — H (lower plot) as 
a function of X t = X® s within Scenario 1, with the same line/color coding as in Fig. [HI 
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Figure 13: Afor </> = h (upper plot) and <f> — H (lower plot) as 
a function of X t = X® s within Scenario 2, with the same line/color coding as in Fig. [HI 
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5 Conclusions 


In this paper we analyzed the scheme dependence of the 0(a t a s ) corrections to the neutral 
CP-even Higgs-boson masses in the MSSM. In a first step we investigated the differences in 
the 0(p 2 a t a s ) corrections as obtained in Refs. [1] and [2]. We have shown that the difference 
can be attributed to different renormalizations of the top-quark mass. In both calculations 
an “on-shell” top-quark mass was employed. The evaluation in Ref. |lj includes the O(e) 
terms of the top-quark mass counterterm, Jmf, however, whereas this contribution was 
omitted in Ref. [2]. We have shown analytically that the terms involving 8m e t do not cancel 
in the 0(p 2 a t a s ) corrections to the renormalized Higgs-boson self-energies (an effect that was 
already observed in the 0(ata s ) corrections in the NMSSM Higgs sector [63]). Numerical 
agreement between Refs. [T] and [2j is found as soon as the 5m e t terms are dropped from 
the calculation in Ref. |Tj. Moreover, as an alternative interpretation, we have shown that 
omitting the 8m £ t terms is equivalent to a redefinition of the finite part of the two-loop field- 
renormalization constant which affects the Higgs-boson mass prediction at the three-loop 
order (apart from a numerically insignificant shift in tan (3 as an input parameter). The 
differences between the two calculations can thus be regarded as an indication of the size of 
the missing momentum-dependent corrections beyond the two-loop level, and reach up to 
several hundred MeV in the case of the light CP-even Higgs-boson. 

In a second step we performed a calculation of the 0(a t a s ) and 0(p 2 a t a s ) corrections 
employing a DR top-quark mass counterterm. We analyzed the numerical difference of the 
Higgs-boson masses evaluated with 8mp s and with 8mf R . By varying the CP-odd Higgs- 
boson mass, Ma , the gluino mass, rrig and the off-diagonal entry in the scalar-top mass 
matrix, X t , we found that in all cases the scheme dependence, in particular of the light 
CP-even Higgs-boson mass, is strongly reduced by going from the full one-loop result to the 
two-loop result including the 0(a t a s ) corrections. The further inclusion of the 0(p 2 a t a s ) 
contributions had a numerically small effect. The differences found at the two-loop level 
indicate that corrections at the three-loop level and beyond, stemming from the top/stop 
sector, are expected at the level of ~ 3 GeV. This is in agreement with existing calculations 
beyond two-loop [371140] . The possibility to use m] )R instead of mf s has been added to 
the FeynHiggs package and allows an improved estimate of the size of missing corrections 
beyond the two-loop order. 
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